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ABSTRACT 

It is pointed out that the total amount of dust in the Universe that is produced in 
steUar evohition in the entire cosmic time is consistent with the observed amount, if 
we add to the dust amount inferred for galactic discs the amount recently uncovered in 
galactic haloes and the surrounding of galaxies in reddening of the quasar light passing 
through the vicinity of galaxies. The inventory concerning the dust closes. This implies 
that dust produced from stars should survive effectively for the cosmic time, and that 
a substantial amount of dust is produced in the burning phase of evolved stars of 
intermedaite mass. 
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1 INTRODUCTION 

A significant amount of gas is ejected into interstellar space 
during the life of stars by stellar wind mass loss and at super- 
nova explosions. Gas contains heavy elements either taken 
from the initial gas or produced during the stellar evolu- 
tion. A significant fraction of heavy elements condenses to 
form dust. Although the mechanisms to produce dust have 
not been well understood, the net amount of dust produced 
from gas with heavy elements can be estimated with some 
confidence. 

Reasonable estimates are available as to how much ma- 
terial is processed by stars by integrating the star formation 
rate as a function of cosmic time. Using recent compilations 
of the star formation rate, we estimate that the total amount 
of material (fuel) processed lies between fifud = 0.004 and 
0.010 (e.g. Fukugita & Peebles 2004, hereafter FP04; Hop- 
kins and Beacom 2006; Fardal et al. 2007; Nagamine et 
al. 2006; Ouchi et al. 2009; Bernardi et al. 2010; Bouwens 
et al. 2010), assuming the Chabrier initial mass function 
(Chabrier 2003). We refer to this amount as the fuel. These 
two numbers represent the curves that pass through close to 
the lower and upper parts of the star formation rate plolQ. 
We remark that the uncertainty seen in the two numbers 
largely arises from the normalisation in low redshift rather 

^ We take as our fiducial the Chabrier initial mass function for 
M < IMq and the Salpeter initial mass function for M > IMq. 
Chabrier takes the mass function slope —2.30 rather than —2.35 of 
Salpeter at a higher mass (he revised to —2.35 in his later version). 
The difference in total mass is about 5%, which is smaller than the 
uncertainty that concerns us here. We also note that the initial 
mass function used in FP04 is close to the Chabrier initial mass 
function, but gives the total mass 10% smaller than the Chabrier. 



than apparently more ambiguous high z behaviours, for in- 
stance, whether the star formation rate declines at z > 3 as 
favoured by recent observations (Ouchi et al. 2009; Bouwens 
et al. 2010). The time span is short at higher redshift and 
the integrated contribution to the total amount of stars from 
2 > 2 is roughly < 20%, so that the high redshift behaviour 
is not important in our argumen10. 

On the other hand, stars and their remnants, white 
dwarfs, neutron stars and black holes, that reside in galaxies 
are estimated to be 



fistar = 0.0030 ± 0.0005 



(1) 



(Fukugita & Peebles 2004; hereafter FP04; 10% upward 
shift is applied to agree with the initial mass function we 
use here). Similar estimates have been made by a number 
of authors (e.g., Shankar et al. 2004; Oohama et al. 2009; 
Bernardi et al. 2010), and they fall in the range indicated 
by the two numbers. There is a gap between the two num- 
bers, Slfuci and Slstar- This may primarily be ascribed to the 
material shed by stars during the evolution, as either stellar 
winds or supernova explosions, whereas there still remains a 
gap between the two values is left as a problem in the future. 

Taking the initial mass function of stars given by 
Chabrier, we estimate that the gas fraction shed by stars 
during the evolution is 0.60 times the mass locked into stars 
and stellar remnants, as summarised in what follows. There- 
fore, the amount of baryons consumed in star formation 
is fifuei = 0.0030 X 1.60 = 0.0048, which agrees with the 

The normalisation at low redshift, which varies among the au- 
thors, is the major source of the uncertainty. Also note that the 



star formation rate increases significantly between 
z = 0.1. 



and 
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amount of fuels estimated from the integration of the star 
formation rate at its lower value of the integral of the star 
formation rate. In the following argument we assume the to- 
tal amount of the fuel is fifuoi = 0.005 for the consistency of 
the argument. This leads to the full consistency of the stellar 
baryon and energy budget, as noted in FP04, including the 
stellar light emission, supernova rates and heavy element 
production, while the accounting must be improved espe- 
cially as to the star formation rate. At the upper edge of the 
star formation rate, the integral of the star formation rate 
seems to overshoot the stars which we see today. 

In our calculation of the dust production we use the 
final mass-initial mass relation for white dwarfs (Serenelli & 
Fukugita 2007; see also Salaris et al. 2009), the mean mass 
of which is O.62A/0 for main sequence stars with the mass 
1 — 8A/0, when averaged over the Salpeter mass function. 
For neutron stars we adopt l.SSM© for the main sequence 
mass 8 — 2^Mq, and 7.5Mq for putative stellar black holes 
for the main sequence mass 25 — lOOM© (Heger et al. 2003). 
The remnant black hole mass is highly uncertain, but the 
fuel that ends with black hole is 8% and the remnant mass 
is only 1.4% the fuel mass with our adopted black hole mass: 
the large uncertainty in the remnant black hole mass is not 
important for our considerations. For completeness let us 
quote that 11% of eq.([T]) are partitioned into white dwarfs, 
8% are substellar and neutron stars are 2%. 

We take this simple model as the base to estimate the 
amount of dust produced in stellar evolution. FP04 have 
calculated that the observed extragalactic background light 
amounts to the energy density Q — 5.1 ± 1.5 x 10~® where 
the optical and far infrared contribute by 2:1, with a 7% 
addition by neutrinos, and that the binding energy in heavy 
element we observe today is altogether Q = — 5.7±1.3x lO"'', 
including nuclei locked in or sequestered from stars. This 
nearly balances the energy in the extragalactic background 
light with the error arising from the two vastly different 
accounting kept in mind. It was also shown that the energy 
output from stars isexpected to be = 5.4 x 10~® for the 
given fuel that amounts to Q, = 0.005. We note that there is 
no missing metal problems at least a,t z — 0, but note that 
80% of heavy elements are locked in white dwarfs and 25% 
of 'metals' is sequestered in neutron stars and black holes. 

In addition, it was also shown (FP04) that 
the present day rate of core collapse supernovae 
0.0079j;!^:|^^3g(100yr Mpc)"^ obtained from the star 
formation rate, integrating over the initial mass func- 
tion from 8Mq to lOOM© agrees with the observed rate 
0.0076l:E!:U[^^o(100yr Mpc)"^ albeit with large errors in 
the observed value. These supernovae produce iron at 
f^Fe = 3.6 X 10~^. A similar estimate of type la supernovae, 
with their normalisation adjusted to the present rate, gives 
the iron abundance 2 x 10~®, which results in the total iron 
abundance flpe = 6 x 10~®, when added, in agreement with 
the estimate of the cosmic iron abundance ftpe = 6.3 x 10~® 
for materials which are not locked up in stellar remnants, or 
sequestered from them. The network of consistencies among 
the numbers shown here points towards the validity of this 
simple framework, concerning the hydrogen fuels consumed 
and the heavy element produced. This tempts us to apply a 
similar consideration to the amount of cosmic dust, which 
was not done in FP04. 

It is yet poorly understood where and how dust is pro- 



duced, and how long does it survive. The calculation we 
show in this paper, based basically only on the final and 
initial mass budgets, circumvents these poorly understood 
aspects and uncertainties. It gives accounting of dust pro- 
duced, which is subject to much less uncertainties. We take 
the Hubble constant Ho = 70km s~^Mpc~^, matter density 
fim = 0.3 in a flat universe whenever necessary. 



2 THE AMOUNT OF COSMIC DUST 

We take the traditional elemental abundance given by 
Grevesse and Sauval (2000), while noting that more recent 
solar abundance estimate by Asplund et al. (2005) leads 
to somewhat a smaller abundance of heavy elements. The 
adoption of the new abundance changes some details of our 
results, but our conclusions are unaffected. The initial so- 
lar abundance is Zq — 0.019 or Z/X\(7, — 0.027 using the 
Grevesse and Sauval table with the present day solar abun- 
dance ^/Xlosurfacc ~ 0.023. with the aid of the solar model 
of Bahcall et al. (2001). The solar abundance would be re- 
duced to Zq — 0.012 if we adopt the revision by Asplund et 
al., 40% smaller than the Grevesse and Sauval value. 

The mean metallicity correlates with luminosity of 
galaxies (Tremonti et al. 2004). When we integrate the 
metallicity over the luminosity function we obtain the cos- 
mic average (Z) = 0.83^©. 

We assume that refractory elements. Si, Fe and Mg in 
interstellar gas or those ejected into interstellar gas all even- 
tually condense to solids. The condensed material is not pre- 
cisely identified, but it is argued that they form effectively 
MgxFe2-xSi04 (the composition of olivine, forsterite or fay- 
alite) with x ~ 1 and MgxFei-xSiOs (enstatite) for a lesser 
amount (Weingartner & Draine 2001). The solar elemen- 
tal abundance of these three refractory elements is similar, 
32—38 ppm per hydrogen, and hence are almost saturated 
if all Si is condensed to the olivine composition. Observa- 
tionally, these elements are known to be highly (> 90%) de- 
pleted from interstellar gas. We assume that these elements 
condense into dust, and they take oxygen by 20% locked up 
in the silicate. 

The other prominent component of dust is carbona- 
ceous material, likely including polycyclic aromatic carbon 
for small grains to graphite for large grains. The estimate for 
the fraction that condense to carbonaceous material varies. 
We take the fraction of condensation to be 50 ±25%, consis- 
tent with the value taken by Weingartner & Draine (2001) 
(Weingartner & Draine 2001). 

Adding the two components we obtain the dust to 
metallicity mass ratio 

r; = dust ~ 0.28 ± 0.07, (2) 

where the error dominantly arises from the two choices of 
the silicate and partly from the assumed amount of carbon 
that condenses into dust. With the solar metallicity for the 
Milky Way, we then obtain 

dust/HI= 1/(135 ±30) (3) 

in agreement with the value adopted by models of dust to 
explain interstellar reddening (Mathis et al. 1977; Weingart- 
ner & Draine 2001) and also with observations. Draine et al. 
(2007) estimated it to be 1/140 for the Milky Way, and find 
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for other galaxies that this fraction varies between 1/100 and 
1/400 with the median 1/190 for SINGS sample of galaxies. 
The fiducial value of 1/100 is usually taken as the dust/HI 
mass ratio in the Milky Way. We assume that the heavy el- 
ements present in the interstellar matter condense into dust 
with the fraction given by eq. ((2|, and dust to HI ratio is 
universal in HI regions. 

The HI survey (Zwaan et al. 2005; Zhang et al. 2008) 
and the H2 survey (Keres et al. 2003) give the cosmic value 
for the sum of atomic and molecular hydrogen to be 

fim+Ha = 5.35 ± 0.87 x 10"" (4) 

where HI : H2 « 0.70 : 0.30. With the assumption that dust 
resides in the HI region and the hydrogen to dust ratio is 
universal we estimate the global dust abundance in galactic 
discs, 

fidust disc = 4.0 ± 1.3 X lO"'' (5) 

This is compared with the dust abundance estimated from 
obscuration due to galaxies in a galaxy survey (Driver et al. 

2007), ^dust galaxy = 3 X 10"^ 

Dust may coagulate to form planets and may be de- 
pleted. The amount we estimated, however, is disturbed lit- 
tle by the planets formation, while they are not entirely 
negligible. Marcy (2005) estimated that 12% of nearby FGK 
stars have detected Jupiter-like planets within 20 AU with 
the planet mass distribution dN/dm ~ . In the recent 
analysis Johnson et al. (2010) indicate that the fraction de- 
pends on the mass of the central star and drops to 3% for M 
dwarfs, which dominate the stars in number. Taking Figure 
4 of Johnson et al. we estimate that 0.052 planets formed per 
1 Mq of fuel consumed. Here we extend the range of stars to 
span the full range of M and A stars, by slightly extending 
the observed range which lies between 0.25 to 2Mq. 

The observed planets are of the Jupiter type, which is 
dominated by hydrogen and helium gas. What concerns us 
here is dust used to form the core of planets. If we take the 
core accretion model for the giant planet formation, the core 
mass is about lOAf® per planet (Mizuno 1980; Pollack et al. 
1996; Rice & Armitage 2003), in agreement with the rocky 
core mass in the solar system planets, which is 8, 10, 12M0 
for Jupiter, Saturn and Uranus (e.g., Lodders & Fegley 1998; 
Guillot 1999). The mass density borne by planetary cores 
that arose from coagulation of dust is fipianet — 8 x 10~^, 
only 1/500 the dust abundance in eq.®. The knowledge 
concerning planets is still immature and the estimates given 
here may be subject to revision in the future. Integrating 
the mass function of planets given above around the Jupiter 
mass, we obtain a rough estimate of the mass density of 
planets to be Jlpianet ~ 7 x 10~^. 

We now calculate the abundance of dust that is ex- 
pected to be produced from the stellar evolution. With the 
Chabrier/Salpeter initial mass function, the gas shed by 
stars is 0.60 times the mass locked into stars with the re- 
sulting remnants given in section 1. The majority of stars 
have solar metallicity. We may consider that this is true even 
at non-zero redshift, possibly except for some galaxies very 
early in structure formation, as indicated by observations 
(e.g., Shapley et al. 2004; de Mello et al. 2004) as well as 
demonstrated by a numerical model for galaxy formation in 
the ACDM universe (Nagamine et al. 2001), and the metal- 
licity estimated from high redshift quasars. The cosmic time 



is short at high redshift where these considerations are more 
relevant. The dominant part of time relevant to star forma- 
tion is at low redshift 2: < 1 — 2. Motivated by these observa- 
tions, we may take the approximation that all stars have the 
normal metallicity the same as at low redshift. We assume 
that the gas shed by stellar evolution has solar metallicity on 
average whether gas taken into star formation is pristine or 
already somewhat enriched. The amount of dust produced 
by gas shed by stars is then 

f^dust produced = ^^fucl X 0.38 X {Z)ri 

= 9.6 X 10"" (6) 

which is about 2.5 times larger than the global amount of 
dust in galactic discs given in ([S]). The total heavy element 
abundance calculated with this approximation correctly re- 
produces the abundance observed at z ~ 0, which in turn 
is constrained by the total energy output as explored by 
extragalactic background light, as shown in FP04. 

We also stress that a number of uncertainties in the 
calculation, such as those in the initial mass function, mean 
metallicity, and the dust/Z ratio, cancel when one compares 
eq. ((5]) with eq. ((6]), so that the error in the relative values 
is not excessively larger. The integral of star formation rate 
still has a significant uncertainty, but it is constrained by the 
observed amount of stars at z = and the fraction of stellar 
mass loss with the aid of the initial final mass relation, as 
well as the energy output argument we quoted in section 1. 
We may take seriously that the gap between the two values, 
expected amount of dust produced and that observed in the 
HI region of galaxies is real. 



3 DUST IN HALOES AND THE CLOSURE OF 
THE DUST ENTRY 

This mismatch between the amount of dust that is ought to 
be produced in stellar evolution and that is observed implies 
that either we miss dust somewhere in the universe or dust 
produced is destroyed and not all survives to now. At the 
beginning, the latter may look natural since lifetime of dust 
is thought to be not very long (Draine & Salpeter 1979; 
Draine 1995; Jones, Tielens & HoUenbach 1996; Dwek 1998). 

We point out, however, that the recent detection of 
dust reddening observed in the galaxy quasar correlation 
promotes us to investigate the former case. Using a large 
quasar sample and yet an even larger galaxy sample of the 
Sloan Digital Sky Survey, taking advantage of its high preci- 
sion multi-colour photometry, Menard et al. (2010, hereafter 
MSFR) have found that the quasar light receives reddening 
when it passes through the vicinity of galaxies. The colour 
dependence of attenuation is in agreement with the dust ex- 
tinction curve known for the Milky Way, although the ratio 
of total to selective extinction Rv ~ 4.9 ±3.2 is not well de- 
termined. This implies that a large amount of dust is present 
around galaxies ranging from 20kpc to a few Mpc, which are 
clearly beyond galactic discs. 

The projected surface density distribution of dust fol- 
lows ~ T-^''-* _ with r the projected distance from the 
centre of the galaxy, similar to the galaxy mass distribution, 
to some 5Mpc. The column density of dust in lines of sight 
is very small but when integrated over a large volume this 
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gives a substantial total amount of dust: MSFR estimated 
the dust mass, 

Mdust 5 X IO^Mq (7) 

for 20kpc < r < rv with ry the virial radius for typical 
galaxies in the sample. 

MSFR inferred fihaio dust ~ 2.1 x 10"'*, but this should 
be taken as a lower limit. For the galaxy sample with me- 
dian redshift of (z) = 0.38 the SDSS observation samples 
galaxies down to ~ 0.25L*. The mean luminosity of the 
sample is then (L) k. Q.7L* . Assuming that the amount of 
dust is proportional to luminosity, the luminosity density 
Lr « 2.2 X W^hLQ (Mpc)"^ means 

^haio dust ~ 2.6 X 10-^ (8) 

if (O is used. MSFR shows that the distribution of dust 
continues to be parallel to that of dark matter associated 
with galaxy to several Mpc, on average beyond halfway 
to neighbouring galaxies. Remembering that the luminos- 
ity density, when multiplied by M/L with M the bound 
mass, gives f» 0.15 and it is compared with « 0.27 of 
global mean mass density, we conclude that 45% of mass is 
present outside the gravitationally-binding radius of galax- 
ies. This is demonstrated using the stacked surface matter 
density distribution derived from weak gravitational lensing 
signal, also in agreement with what we expect for the dark 
matter distribution in the CDM universe (Masaki, Fukugita 
& Yoshida 2011, in preparation). This suggestss that the 
amount of dust we estimated must be multiplied by 1.8, in- 
cluding the distribution beyond the virial radius, to obtain 
the global abundance, if dust follows the dark matter dis- 
tribution as observationally indicated in MSFR. This means 
that figaiaxy+vicinity dust ~ 4.7 X lO"'', if the viciuity of galax- 
ies beyond the gravitationally bound region is included. 

The uncertainty in the estimate of the amount of dust 
outside the galaxies is admittedly large, but the order of the 
amount we see indicates that the amount of dust in galactic 
discs must be almost doubled for the total amount. This 
means that the total amount of dust 

^dust ~ 9 X 10-^ (9) 

which is close to what is expected from stellar evolution and 
the amount of cooked fuels Q. 

It has been discussed that dust is destroyed efficiently 
by sputtering in the hot environment with the time scale, 
t « 10^yr(nH/cm~^)"^(a/0.1^im) for T^g ~ 10®, the virial 
temperature of L* galaxies (Draine & Salpeter 1979). In the 
halo environment, nu ~ 2 x 10~* at lOOkpc (Fukugita & 
Peebles 2006), so that lifetime is longer than the age of the 
universe. It seems likely that dust in galactic haloes trans- 
ported there from the galaxy during the course of its forma- 
tion and evolution may survive for the cosmological time. 

Our argument indicates that the longevity of dust 
should also apply to galactic discs. The total amount of dust 
we obtained is what is ought to be produced. We have no 
more fuels to yield extra dust. If dust would be destroyed for 
some reasons, it must be replenished by regeneration from 
interstellar matter. 

We also remark that dust may be efficiently transported 
to outside galaxies by galactic winds. For some galaxy sam- 
ples it is observed that a substantial fraction of the mass of 
formed stars is outflowed by galactic winds (e.g., Heckman 



et al. 2002; Pettini et al. 2002; Veilleux et al. 2005; Rupke 
et al. 2005; Weiner et al. 2009), as also supported by plau- 
sible arguments and simulations (e.g., Madau et al. 2001; 
Aguirre et al. 2001; Zu et al.). Whether a large part of dust 
is present in some clumps, such as absorbing clouds (Menard 
et al, 2007), raises an interesting problem, but in the present 
study, where we are concerned with the coarse-grained mean 
abundance, we have no resolution to this problem. 



4 CONCLUSIONS AND DISCUSSION 

We have shown that the amount of dust expected in stellar 
evolution in the cosmic time agrees with what we observe, 
when we take account of dust that is uncovered in galac- 
tic haloes and in the vicinity of galaxies from reddening of 
quasar light passing through nearby galaxies (MSFR), which 
roughly doubles the total amount when added to that ob- 
served in galactic discs. Namely, we are observing all dust 
produced in stellar evolution, which means that dominant 
portion of dust should survive, at least effectively for cosmic 
time. This contrasts to the conventional thought that life- 
time of dust is much shorter than the cosmic time. This does 
not imply, however, that dust be intact. If it is destroyed in 
some processes in interstellar or intergalactic space, it must 
be regenerated from the gas. We emphasise that we have no 
extra fuels to replenish destroyed dust, meaning the closure 
of the cosmic energy inventory, which in turn places a sig- 
nificant constraint on the physics of dust. The dust amount 
used to form cores of planets are only 1/500 in the disc. 

If the dust production were ascribed solely to core col- 
lapse supernovae, the total amount of dust we estimated 
means the productin of O.2M0 per core collapse. This mass 
is consistent with the prediction in some theoretical calcula- 
tions (e.g., Kozasa et al. 2009; Todini & Ferrara 2001), but 
is 2 orders of magnitudes larger than is actually observed 
at type II supernovae, 10~* — 10~"^Mq. In our considera- 
tion we do not necessarily mean that dust is produced at 
the time or just after supernovae. It may in part arise from 
mass loss in the giant star phase before supernovae, or from 
gas ejected by supernovae that may condense to dust later 
in the interstellar space (Draine 2009). 

The closure of the inventory tempts us to look at the 
partition more closely as to the origin of dust. From the con- 
sideration of the final versus initial stellar masses we expect 
that the stellar wind of stars with 1 — 8Mq produces 56 % of 
dust, where the dominant part arises from the AGB phase 
(see Serenelli & Fukugita 2007). Stars that end with core 
collapse supernovae produce 43 % of dust, including both 
supernovae (and after supernovae) and mass loss in the pre- 
supernova phase. This means that the amount of dust pro- 
duced in core collapse is at most O.OSMq, smaller than the 
theoretical estimate cited above. Normalising the type la su- 
pernova rate to the observed value at redshift z « as in 
FP04, we infer that type la supernova disrupt 5% of white 
dwarfs which contribute 1% of dust at the time or some time 
after the explosion. There is no evidence observed that type 
la supernovae yield dust. 

Our consideration suggests that a substantial fraction 
of dust is generated from the material lost in the burning 
phase of evolved stars. The inventory gives useful circum- 
stantial constraint on the production and evolution of dust. 
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although it does not directly suggest anything concerning 
the mechanisms. 
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